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We derive a general scaling relation for the anomalous Hall effect in ferromagnetic metals involving 
multiple competing scattering mechanisms, described by a quadratic hypersurface in the space 
spanned by the partial resistivities. We also present experimental Hndings, which show strong 
deviation from previously found scaling forms when different scattering mechanism compete in 
strength but can be nicely explained by our theory. 


Despite recent success in recognizing the importance 
of intrinsic Berry curvature in the anomalous Hall ef¬ 
fect (AHE) [1-7], there remains a basic conflict between 
theories and experiments on the scaling behavior of the 
anomalous Hall resistivity [8]. Working under the as¬ 
sumption of a single type of scatterers, theories invari¬ 
ably show that the anomalous Hall resistivity is a sum 
of linear and quadratic terms in the longitudinal resistiv¬ 
ity, with the former identified as from skew scattering [9] 
and with the latter as from both intrinsic and side jump 
mechanisms [10]. However, many experiments do not fol¬ 
low such a simple scaling relationship. 

Recent experimental [7, 11, 12] and theoretical stud¬ 
ies [13, 14] indicate that a major cause of the breakdown 
of the simple scaling is the presence of more than one 
source of scatterers. By working with molecular beam 
epitaxial films, Tian et al. [7] were able to limit scatter¬ 
ing of electrons to two sources, one by interface roughness 
and one by phonons, with independent control on their 
strengths through the film thickness and sample temper¬ 
ature. By taking into account the distinction of scaling 
parameters with respect to different scatterings, a new 
phenomenological scaling form were found which fits the 
experimental data much better. However, its theoretical 
derivation has been lacking so far. 

In this Letter, we derive the general scaling form of the 
anomalous Hall resistivity as a function of multiple com¬ 
peting scattering mechanisms. This is also motivated by 
the new experimental data to be presented in this work, 
which deviates from our earlier results in the regimes 
when two sources of scatterings compete in strength. We 
find that the anomalous Hall resistivity is generally a 
quadratic surface passing through origin in the space of 
partial longitudinal resistivities. This reduces to the sim¬ 
ple one-variable scaling along each axis, but contains an 
extra parameter that describes competition between the 
two scatterings. 

The general scaling relation obtained from our theory 
with multiple competing scatterings can be written in the 


following form 

Pah = (^Pxx Y ^ ' CiPiPxx + ^ ^ CijPiPj + ^ ^ ctiPi, (1) 

i ij iGS 

where the Hall resistivity pah is expressed using the 
partial longitudinal resistivities pi^) as scaling variables 
which are changed through the disorder concentrations 
of each scattering source, i and j label different scatter¬ 
ing sources, the coefficient c encodes the intrinsic Berry 
curvature contribution, and coefficients c;, Cij, and a; re¬ 
main constant when disorder concentrations are changed 
(they do depend on the source of scattering as indicated 
by the indices), and i G S indicates that only static dis¬ 
order scatterings are included in the summation. With 
Matthiessen’s rule pxx = ^ can be seen that ge¬ 

ometrically PAH spans a quadratic hypersurface in the 
space of Pi’s. Along each p^-axis, the scaling takes the 
form of a parabolic curve a^pi -|- (c -I- -|- Cii)p^ passing 

through the origin, which is reminiscent of the old simple 
scaling [8] (which applies only for a single type of scat¬ 
tering) with the linear term (a;) from skew scattering 
and the quadratic term from both intrinsic (c) and side 
jump (ci + Cii) mechanisms. Furthermore, cross terms of 
PiPj generally exist, manifesting the competition between 
different scatterings. 

Sketch of theoretical derivation .— In the following, we 
provide a justification of our general scaling formula from 
microscopic theory. Although the different mechanisms 
have been envisaged from the early theories of AHE, the¬ 
ories which could systematically account for all the mech¬ 
anisms were developed relatively recently [8] . It has been 
shown that the different approaches are equivalent in the 
good metal regime with ep-r ^ 1 (e^? is the Fermi energy 
and T is the scattering time) which is what we are are 
concerned with here. 

Theoretically, it is more convenient to calculate con¬ 
ductivity which, in Kubo formalism, is related to the 
equilibrium current-current correlation function and can 
be casted into a familiar form of the Kubo-Streda for¬ 
mula [15]. There are two terms for the Hall conduc- 
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FIG. 1. (a) Correction of velocity vertex by a ladder diagram. 
The dressed vertex v contains factor jr. (b-d) correspond 
to our defined groups 1 to 3 respectively. The colors of two 
Green function lines at each vertex indicate whether they be¬ 
long to the same band. The dependence of leading order terms 
of each part on longitudinal resistivities is indicated. 

tivity (TAH, one of which (ct^h) i® states near 

the Fermi surface which contains all the important scat¬ 
tering effects, and the other term crJJjj involves all the 
occupied states below the Fermi level and only con¬ 
tributes to the intrinsic part. Hence to assess the 
scattering effects, one only needs to focus on = 

2 ^Tr(Da;G^(£:F)wyG^(eF))c- Here A is the system size, 
Vx and Vy are the velocity operators, and are the 
retarded and advanced Green functions respectively, and 
indicates a disorder configuration average. The 
usual practice is to treat scattering as perturbation and 
do expansion in the small parameter For sim¬ 

ple models with single type of scattering, one can identify 
relevant diagrams corresponding to intrinsic, side jump, 
and skew scattering contributions, and group the terms 
according to their dependence on disorder density rii, as 
discussed in detail by Sinitsyn et al. [16] and reviewed 
by Nagaosa et al. [8]. Because both skew and side jump 
contributions depend on the detailed form or composition 
of scattering potentials, the variation of which along the 
change of experimental parameters can naturally cause 
the complication, especially when more than one source 
of scatterings are competing. However, since our focus is 
on the scaling relation, it is sufficient to extract the gen¬ 
eral structure and to link ctah to the longitudinal trans¬ 
port coefficient (which in Kubo formalism takes similar 
form as by replacing Vy with Vx) without perform¬ 
ing detailed calculations. Our analysis below applies to 
general multiband systems with multiple sources of scat¬ 
terings. 

For longitudinal transport in metals, it is well-known 
that the scattering has two important effects [17]: first, 
it introduces a finite lifetime for quasiparticles through 
self-energy correction to the Green functions; second, it 
corrects the velocity vertex, making axx cx where 


T*'' is the transport scattering time. In the so-called lad¬ 
der approximation, the velocity vertex is dressed by a 
ladder diagram as shown in Fig. 2(a). Note that only in¬ 
traband velocity vertex needs to be dressed, because the 
product involved, [G^\nk[G^]mki is large (ex r) to com¬ 
pensate the small factor from adding a scattering line 
(oc I/t) only when the band indices m and n are equal; 
while for n ^ m, the product can be at most regarded 
as a factor of r°. The renormalized velocity vertex v 
acquires a factor of jr, and when combined with the 
two Green function legs, one has [G^]nkVnn[G^]nk oc r*'' 
upon summation over states \nk). When multiple sources 
of scatterings are present, r*'' (r) would correspond to 
the total transport (usual) scattering time. Assuming 
there is no correlation between different sources of scat¬ 
terings, then each scattering contributes to t*"' (and r) 
independently as dictated by Matthiessen’s rule, i.e., 
1/r*'' = where i labels the source of scattering, 

and hence pxx =J2iPi- 

One notes that the same structure appears in 
as well. However, unlike axx, the interband coherence 
effects including the interband velocity component and 
the interband off-shell scattering process (i.e. virtual 
transitions to states off the Fermi surface) are impor¬ 
tant [8], such that in the leading order diagrams each 
velocity vertex can be either intraband (dressed) or in¬ 
terband (not dressed). This suggests us to group the 
terms in tr^jj according to the characters of velocity ver¬ 
tices involved, which naturally leads to three groups (see 
Fig.2(b-d)) [19], i.e., diagrams with 1) two interband ver¬ 
tices; 2) one interband vertex and one intraband vertex; 
and 3) two intraband vertices. Within each group, we 
can then identify and analyze the leading order terms. 

Following this scheme, for group 1, the leading term 
is the single bubble without any off-shell scattering lines. 
This combined with crj.^ gives the intrinsic contribution. 
For group 2, the part of dressed intraband velocity ver¬ 
tex along with two Green function legs, according to our 
previous discussion, is oc r*'' cx axx ~ Pxx in terms of 
its scattering dependence. Meanwhile, its counterpart at 
the other vertex is of zeroth order in t. For the remaining 
middle part as indicated in the dashed box in Fig. 2(c), 
there must be at least one off-shell scattering line which 
brings a small factor oc l/r, oc pi for a scattering of type 
i. The important point is that this group of terms share 
a common factor IjPxx from the vertex part. Hence sum¬ 
ming over i leads to a generic form of 'Ylii^iPilPxx with 
coefficients Ci that do not depend on the disorder con¬ 
centration. Similar analysis applies to group 3. The two 
vertex parts each brings a factor of IjPxx- In the middle 
part there must be at least two off-shell scattering lines 
(~ PiPj) or one three-scattering line (~ pi) [20]. Note 
that the latter does not occur for dynamic disorder scat¬ 
tering (such as phonon) [14, 18, 21] and it corresponds 
to the conventional skew scattering. 

Collecting all the leading order terms and expressing 
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them in terms of the longitudinal partial resistivities, we 
obtain the result 

C + 'y ' CiPijPxx “b Cy ' C-ijPiPj + y ' UiPi) jPxx- 
i ij ieS 

( 2 ) 

The c term here corresponds to the usual intrinsic con¬ 
tribution. Transferring to the Hall resistivity pah — 
—tTAHPxa: [ 22 ], we thus obtain the general scaling rela¬ 
tion in Eq.(l). 

According to our experimental setup, if we specialize 
Eq.(l) to two major scattering mechanisms: one static 
and one dynamic, which could be tuned experimentally 
into competition, then we have pah = cp^^: + {copxxo + 
^iPxxT)Pxx T (^^OOPxxQ ~b ^QIPxxQPxxT “b ^llPxxT ~b ^Pxxo'ji 
where pxxo is the residual resistivity due to static impu¬ 
rities at low temperatures and Pxxt{= Pxx — Pxxo) is due 
to dynamic disorders (mainly phonons at higher temper¬ 
atures). Since there are only two independent variables 
out of the three {pxx,PxxO, Pxxt), if we choose to write in 
term of the two partial resistivities PxxO and PxxT, 

Pah = oipxxo + PoPxxo + IPxxoPxxT + PiPxxT- ( 3 ) 

with /3o = c + Co + Coo, 7 = 2c -I- co -b ci -I- coi, and 
/3i = c -\- Cl + Cii- One notes that there are only four 
scaling parameters: a, /?o, /3i, and 7 , which can be 
extracted from experiment. Except for the skew scat¬ 
tering term apxxo (which is absent for dynamic disor¬ 
ders [7, 14, 18, 21]), the scaling is bilinear in the two 
partial resistivities, and can be viewed as a nice combi¬ 
nation of simple scalings of each type of scattering plus 
a cross term due to competitions. Furthermore the coef¬ 
ficients j3o and /?i of the two quadratic terms contain not 
only the intrinsic Berry curvature contribution c but also 
side-jump contributions due to each scattering. In the 
following, we shall see that this result not only justifies 
the previous phenomenological scaling in the low con¬ 
ductivity regime [7], but also nicely captures its strong 
deviations in the high conductivity regime. 

Experimental evidence. — Now we proceed to the dis¬ 
cussion of our experimental data. Only after the previous 
work of Fe on GaAs(OOl) [7], we find that the epitaxial 
Fe films on MgO(OOl) have much higher quality, as re¬ 
flected by a much higher electrical conductivity for a film 
at the same thickness and temperature. It is from these 
high quality Fe film data that we discover significant de¬ 
viations from the scaling proposed in [7] , which call for a 
more rigorous theory for the scaling of ABE. In Fig. 2(a) 
and Fig. 2(b), the longitudinal resistivity and anomalous 
Hall resistivity measured from 5 K to 310 K are shown re¬ 
spectively, from which we can obtain the anomalous Hall 
conductivity ctah and plot it against a'^x ia Fig. 2(c). Ac¬ 
cording to the scaling ctah = {aa~xo + li(^xxo)(^lx + & 
obtained phenomenologically in [7] (hereafter referred to 
as the TYJ scaling), a linear dependence should have 
been expected, however, except for the 8 nm data, all 


the curves here show obvious deviation from the linear 
relation. Nevertheless, once zooming in the data in the 
low conductivity regime (where most of the data from Fe 
on GaAs located, see. Fig.4 in Ref. [7]) as given in the 
inset, we can clearly see that they are indeed well de¬ 
scribed by the TYJ scaling. Gomparing the present data 
with that of Fe on GaAs(OOl), we find that the break¬ 
down of the TYJ scaling occurs in the high conductivity 
regime while in the real clean limit only the skew scat¬ 
tering term counts [8]. In fact, by a closer look one could 
also find its traces in Fig.4 of Ref. [7], although not so 
obvious as here in Fig. 2(c). 

Realizing the limitation of the scaling in Ref. [7], we 
now test whether our general scaling relation works here. 
For better comparison and analysis, rewriting Eq.(3) us¬ 
ing conductivities Uxxo and Oxx gives 

- ctah = oia^^^alx + (/^o + /3i - l){'^xxo^xx)'^ 

+ (7 - 2/3i)(T,j,jo(T2,a; -I- /3i. (4) 

Gompared with the TYJ scaling, one observes that the 
third term on the right hand side of Eq.(4) is formally 
new. We will expose the deeper connection between the 
two in a while. To extract the four fitting parameters, 
first, by using the experimental data at 5 K in Fig. 2, 
we plot the residual conductivity ctaho versus cTxxO for 
different film thicknesses in Fig. 2(d). Because of the 
limited thermal excitations at 5 K, Gxx — ctxxOj Eq.(4) 
becomes —ctaho = ckcTxxO + /3o for different film thick¬ 
nesses. Fitting the data using this relation (see Fig. 2(d)), 
we could first pin down the skew scattering coefficient 
a = -1.49 X 10-3. 

Then, by using the experimental data with variable 
temperature from 5 K to 300 K in Fig. 2, we show in 
Fig. 2(e) a series of —uah vs Oxx plot for different film 
thicknesses. Here we subtract the skew scattering contri¬ 
bution from —ctah in Eq.(4), and use it to fit the corre¬ 
sponding experimental data, as shown by the solid lines 
in different colors in Fig. 2(e). On the other hand, one 
notes that after the subtraction of the skew term, the 
right hand side of Eq.(4) is in fact a parabolic function 
of (Jxxo^xx- Indeed, if plotting (—ctah — ota^xo^xx) ver¬ 
sus fyxxCTxxo nsing the experimental raw data for different 
film thicknesses in Fig. 2(f), we find strikingly that all the 
curves from Fig. 2(e) collapse on top of each other almost 
perfectly as indicated by the eye-guided line. 

Plotting in Fig. 2(g) all the remaining three fitting pa¬ 
rameters obtained in Fig. 2(e), we find that the three pa¬ 
rameters are of the same order of magnitude and their 
values are quite stable across different film thicknesses. 
Particularly, j3o and /3i (which include intrinsic and side 
jump) are nearly constant. Here we have two sepa¬ 
rate ways to extract Po'- either from the interceptions of 
Fig. 2(d) using the experimental data at 5K or from the 
fitting of each curve in Fig. 2(e) using all the data from 
5 K to 300 K. As shown in Fig. 2(g), the obtained val- 
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FIG. 2. (a,b) Temperature dependence of pxx and pah of epitaxial Fe films on MgO(OOl), respectively, (c) ctah versus 
for Fe films with various thicknesses. The inset shows the regime in which the data can be well fitted with the TYJ scaling, 
while systematic deviations from the TYJ scaling develop in the larger range, (d) ctaho versus at 5 K. The red line 
is the linear ht and the slope a is the skew scattering coefhcent. (e) Anomalous Hall conductivity with the skew scattering 
contribution subtracted (— ctah — CKa~^Qa^^) versus a^x for various film thicknesses. Fitting curves are based Eq.(4). (f) 
(—UAH — aa~^gaxx) versus o-xxO'~^g for different film thicknesses which almost all collapse onto a single curve. The red line is 
a guide of eye. (g) Fitting results for parameters /3o, /3i, and 7 . The difference of the values of /Jq obtained from the fitting in 
present figure (brown dots) and from the interception in (d) is within the error bar. 



FIG. 3. The quadratic surface of pah(PxxO, Pxxt) according 
to Eq.(3) with the values of four parameters determined in 
Fig. 2. The deviation of each raw data point from the surface 
is within the size of the dots. 

ues of /3o agree with each other very well. 7 has a slight 
deviation for small thicknesses which could be originated 
from the finite size effect of dynamical scatterings when 
the film is very thin. Importantly, all three parameters 
quickly approach constant (bulk) values when the film 
thickness is greater than 20 nm. 

Now with all four parameters a, /3o, /3i, and 7 
fitted from experiment, we plot in Fig. 3 the surface 


Payi{PxxO,Pxxt) as described by our general scaling 
Eq.(3). Remarkably, one observes that almost all the 
data points lie perfectly on this paraboloid surface, con¬ 
firming the validity of our scaling relation. 

Discussions .— We point out the important connection 
and contrast between our general scaling and the TYJ 
scaling in [7]. We rearrange Eq.(4) into the following 
form 

- ctah = -I- (/3o - Pi)cr^^Qal^ + 

( 7 - 2 / 3 i )( '^xxO^xx ^xxO^xx )• ( 5 ) 

While the first three terms on the right hand side have 
the same form as the TYJ scaling, one notes that the 
last term is negligible in both the high temperature 
limit (where cr^xo 3> (Txx) and the low temperature limit 
(where axx — cTxxo)- This hence explains the success of 
the phenomenological scaling as for Fe films on GaAs [7] 
as well as in the inset of Fig. 2(c), due to the smallness of 
the last term in this regime. Accordingly, for each film 
thickness in Fig. 2(c), by connecting with dash line the 
two data points at the highest and lowest temperatures 
respectively, we are in fact describing the correspond¬ 
ing AHE using the first three terms of Eq.(5). Obvi¬ 
ously it significantly deviates from the reality where the 
real experimental points are connected by the eye-guided 
line, especially for thicker films with larger conductiv¬ 
ity or relatively weaker static impurity scattring. There¬ 
fore, it is in the intermediate temperature regime where 
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the static and dynamic impurity scatterings are compet¬ 
ing with each other for the AHE, that the last term of 
Eq.(5) manifests itself unambiguously as a critical fac¬ 
tor to the AHE (one notes that the coefficient of the last 
term (7 —2/3i) = cq —ci-fcoi—2cii is purely of scattering 
effect and is consistent with its interpretation as from the 
competition between scatterings). In view of the new ex¬ 
perimental facts and deeper understanding presented in 
this work, the ‘proper scaling’ found phenomenologically 
in the past [7] is no longer appropriate. 

Finally it should be noted that the coefficients of 
the second and the third term of Eq.(5): /3o — /3i = 
Co — Cl -|- Coo ~ Cii and Pi = c + ci -\- cn also have new 
meaning compared with the scaling in [7]. One observes 
that (/3o — Pi) is also purely of scattering effects and the 
only presence of the intrinsic Berry curvature contribu¬ 
tion c in Eq.(5) is in Pi. Our work demonstrates that 
with independent control of scatterings, it is possible to 
experimentally extract ai and (c + Ci + cu) for each scat¬ 
tering source. Nevertheless, as in both Eq.(5) and Eq.(3), 
the intrinsic contribution c is entangled with the side- 
jump contributions (c^ -|- cu). Given that only the four 
parameters (and their linear combinations) can be ex¬ 
tracted, therefore it seems challenging to unambiguously 
separate out the intrinsic part c from such experimental 
approach. This will surely motivate further experimental 
studies in near future. 
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